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Introduction 

Cyclic adenosine 3', 5 '-monophosphate (cAMP) is a 
widely used biochemical messenger, transducing extra- 
cellular stimuli into a myriad of cellular responses. In the 
late 1970s, data emerged showing that hormone receptors 
acting via cAMP have distinct effects on cardiac function. 
At hormone levels where (3-adrenergic and prostaglan- 
din receptors induced similar changes in overall cAMP 
and PKA activity, only (3-adrenergic receptors were found 
to increase left ventricular contractility, glycogen me- 
tabolism, and troponin I phosphorylation (Brunton et al., 
1979; Hayes et al., 1979) . This indicated that cAMP does 
not indiscriminately activate its downstream targets, but 
rather it acts in a context-dependent manner. 

To explain the context-dependent specificity of cAMP, 
Brunton et al. (1981) proposed the compartmentation 
hypothesis, in which receptor-specific cAMP micro- 
domains are distinctly coupled to cellular functions. 
Although this concept was initially quite controversial 
(Steinberg and Brunton, 2001), a variety of technologies 
has been developed over the past two decades that have 
enabled direct evidence of cAMP compartmentation. A 
key focus has become to understand the biochemical 
and biophysical mechanisms underlying cAMP com- 
partmentation. Because of the limited ability to specifi- 
cally perturb and measure all aspects of cAMP signaling 
experimentally, computational models have been devel- 
oped to help understand how cAMP compartments are 
formed. Indeed, computational models are well suited 
for identifying biological mechanisms, predicting down- 
stream consequences, and reducing the complexity of 
large datasets (Yang and Saucerman, 2011). 

As the experimental efforts to measure and manipulate 
cAMP compartmentation have been well reviewed else- 
where (Steinberg and Brunton, 2001; Saucerman and 
McCulloch, 2006; Willoughby and Cooper, 2007; Karpen, 
2014; Rich et al, 2014) , this Perspective will concentrate 
on the specific insights into cAMP compartmentation 
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provided by computational models. Computational mod- 
els have been used to evaluate a range of potential cAMP 
compartmentation mechanisms: localized cAMP syn- 
thesis, localized cAMP degradation, physical barriers to 
diffusion, cAMP buffering, cell shape, and cAMP export 
(see Fig. 1). After briefly summarizing key motivating 
experimental measurements, we will describe model 
predictions related to each of these potential mecha- 
nisms. We will then discuss future directions including 
necessary experimental validations of key model predic- 
tions and the incorporation of cAMP compartmenta- 
tion into multi-scale computational models. 

Experimental measurements of cAMP compartmentation 
Biochemical approaches. The initial measurements of 
cAMP compartmentation were performed by cellular frac- 
tionation and radioimmunoassay. Corbin et al. (1977) iso- 
lated particulate and soluble fractions of rabbit heart 
homogenates, finding that about half of the total cAMP 
content was bound to PKA regulatory subunit in the par- 
ticulate fraction. Increasing cAMP synthesis or blocking 
its degradation caused disproportionate [cAMP] increases 
in the soluble fraction (Corbin et al., 1977). Although 
activation of both (3-adrenergic and prostaglandin re- 
ceptors increased soluble cAMP and PKA activity in 
heart homogenates, only (3-adrenergic receptors ele- 
vated cAMP and PKA in the particulate fraction (Hayes 
et al., 1980) and triggered downstream increases in con- 
tractility and glycogen metabolism (Brunton et al., 1979) . 
A limitation to these biochemical approaches is that they 
destroy the intact cellular environment, and particulate 
fractions contain a wide range of membranes, sarcomeres, 
and organelles. 

Electrophysiological approaches. Creative use of patch- 
clamp electrophysiology allowed more direct measure- 
ment of cAMP compartmentation in live cells. Jurevicius 
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and Fischmeister (1996) used a microperfusion system, 
finding that local application of the adenylyl cyclase 
(AC) agonist forskolin enhanced L-type Ca 2+ currents 
globally, whereas locally applied |3-adrenergic agonist 
isoproterenol produced only local elevations in L-type 
Ca 2+ currents. These approaches were further enhanced 
by the use of CNG channels. Rich et al. (2000) used 
patch clamp of HEK293 cells expressing cAMP-sensitive 
CNG channels, finding that forskolin induced much 
higher submembrane [cAMP] than global [cAMP] . 

Fluorescent biosensors. A wide range of fluorescent biosen- 
sors for cAMP has been engineered. The first used fluo- 
rescein and rhodamine-labeled regulatory and catalytic 
subunits of PKA, where cAMP binding lead to a decrease 
in fluorescence resonance energy transfer between the 
fluorophores, allowing visualization of [cAMP] gradients 
induced by serotonin (Bacskai etal., 1993). Zaccolo etal. 
(2000) improved on this approach by fusing regulatory 
and catalytic subunits of PKA with cyan and yellow fluo- 
rescent proteins, creating a genetically encoded PKA- 
based biosensor. Their biosensor was used to visualize 
micrometer-scale cAMP gradients induced by (3-adrener- 
gic agonist in live cardiac myocytes (Zaccolo and Pozzan, 
2002) . Alternative cAMP biosensors have used conforma- 
tional changes in the cAMP-binding protein Epac (DiPilato 
et al., 2004; Nikolaev et al., 2004) or the cAMP-binding 
domain of the hyperpolarization-activated CNG channel 
2 (termed HCN2-camps) (Nikolaev etal., 2006). 

Localized cAMP degradation 

By far the most prominently recognized mechanism for 
cAMP compartmentation is localized cAMP degradation 



by phosphodiesterases (PDEs) (Fig. 1 A) (Francis et al., 
2011). Jurevicius and Fischmeister (1996) provided the 
first evidence of PDE-mediated cAMP compartmentation, 
showing that PDE inhibition allowed local (3-adrenergic 
stimulation to enhance Ca 2+ currents globally in frog 
ventricular myocytes. Inhibition of PDEs ablated com- 
partment-specific cAMP dynamics and receptor-medi- 
ated cAMP gradients when using cAMP-sensitive CNG 
channels (Rich et al., 2001) or PKA-based fluorescent 
biosensors (Zaccolo and Pozzan, 2002). 

Nearly all computational models have predicted a 
quantitatively important role of PDE-mediated cAMP 
degradation in the formation of cAMP gradients (see 
Table 1). The model developed by Rich et al. (2001) 
predicted that prostaglandin El (PGEl)-stimulated PDE 
activity was critical for quantitatively explaining the 
transient plasma membrane cAMP signals measured ex- 
perimentally. A subsequent experimental and compu- 
tational study from their group further demonstrated 
that the enhanced PDE activity was caused by PKA-me- 
diated PDE phosphorylation (Rich et al., 2007) , agreeing 
with previous experimental evidence of negative feed- 
back by PKA on [cAMP] (Rochais et al., 2004) . Saucerman 
et al. (2006) modeled cAMP and PKA signaling with im- 
age-based cardiac myocyte geometry, predicting that 
PDE regulated the magnitude of cytosolic cAMP/PKA gra- 
dients. Simulations with realistic neuronal geometries 
(Neves et al., 2008) have also explained how PDEs con- 
tribute to experimentally observed cAMP gradients be- 
tween dendrites and the cell body (Bacskai et al., 1993) . 

The expression of multiple PDE isoforms and their 
localization to various cellular structures suggest that 
PDEs may regulate distinct cAMP compartments. In 
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Figure 1 . Predicted mecha- 
nisms of cAMP compartmen- 
tation. (A) PDEs can locally 
degrade cAMP to create gra- 
dients. (B) cAMP synthesis by 
AC can elevate local [cAMP] . 
(C) Physical barriers restrict 
cAMP diffusion. (D) cAMP 
binding to PKA can reduce 
the freely diffusing [cAMP]. 

(E) Cell shapes that alter the 
surface-to-volume ratio can 
alter the local balance of cAMP 
synthesis and degradation. 

(F) Export of cAMP from the 
cell by MRPs can decrease 
local [cAMP]. 
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most cells, PDE3 and PDE4 are the dominant PDE iso- 
forms for cAMP degradation (Francis et al., 2011). Al- 
though PDE3 contains membrane-binding domains and 
appears to target to intracellular membranes, PDE4 is 
primarily bound to A kinase-anchoring proteins (AKAPs) 
(Francis et al., 2011; Kapiloff et al., 2014) and localized 
in striated patterns, colocalized with either sarcomeric 
M-line or Z-lines (Mongillo et al., 2004). Measurement 
of cAMP signals using PKA-based fluorescence resonance 
energy transfer biosensors revealed that PDE4 dominates 
regulation of both basal and (3-adrenergic-stimulated 
cAMP, whereas PDE3 dominated cAMP responses to 
forskolin in cardiac myocytes (Mongillo et al., 2004) . 
Similarly, PDE4 plays a dominant role regulating sub- 
sarcolemmal cAMP, as measured by CNG channels ex- 
pressed in cardiac myocytes (Rochais et al., 2004). 

Several models have examined how PDE localization 
contributes to cAMP compartmentation. Saucerman et al. 
(2006) measured and modeled (3-adrenergic-stimulated 
gradients of cAMP/PKA signaling between sarcolemma 
and cytosol, predicting that the measured time delays 
from sarcolemmal to cytosolic PKA activity were insensi- 
tive to PDE localization or mobility unless cAMP diffu- 
sion rates were slowed. Iancu et al. (2007) developed a 
multi-compartmental model of cAMP signaling in car- 
diac myocytes, including PDE2, PDE3, and PDE4 iso- 
forms, and distinct sarcolemmal domains. Their model 
predicted that heterogeneous PDE localization (not the 
PDE isoform per se) was critical for maintaining cAMP 
gradients under basal and stimulated conditions (Iancu 
et al., 2007). This is a helpful model prediction because 
it is difficult to experimentally dissect PDE isoform 
localization from their distinct kinetic and regulatory 
properties. Heijman etal. (2011) further extended this 
model (as discussed below) , predicting that PDE3 plays 
the dominant role in regulation of baseline caveolar 
and cytosolic [cAMP]. Chen et al. (2008) developed 
models with varying localization of PDEs, showing that 



colocalization of PDEs with ACs provided optimal com- 
partmentation. In contrast, they predicted that if PDEs 
were uniformly distributed, PDEs expressed at bio- 
chemically reasonable values would be insufficient to 
generate cAMP gradients (Chen etal., 2008). 

Although most models of cAMP compartmentation 
have been deterministic (repetition of the same simula- 
tion gives the same results), Oliveira et al. (2010) devel- 
oped a stochastic cAMP model to reflect the variability 
caused by small numbers of signaling molecules. Using 
this model, they predicted that PKA-mediated phos- 
phorylation of cytosolic PDE4D creates a negative feed- 
back critical for formation of submembrane cAMP 
microdomains in HEK293 cells (Oliveira et al., 2010), 
consistent with prior experiments (Rich et al., 2001; 
Terrin et al., 2006). 

In summary, a wide range of experiments and models 
supports the key role of PDEs as necessary for generation 
of [cAMP] gradients. However, some of these models 
differ considerably in the extent to which PDE-mediated 
cAMP degradation is sufficient for [cAMP] gradients, 
or whether other mechanisms are required as well. 

Localized cAMP synthesis 

cAMP is synthesized by ACs from ATP. There are several 
types of ACs that differ in their distinct patterns of molecu- 
lar regulation and tissue-specific expression (Hanoune 
and Defer, 2001 ; Willoughby and Cooper, 2007) . Although 
all transmembrane ACs are activated by the binding of 
GTP-bound G sa or the AC agonist forskolin, they differ 
in their molecular regulation and expression (Hanoune 
and Defer, 2001; Willoughby and Cooper, 2007). Solu- 
ble ACs have also been identified, which are activated 
by bicarbonate, and several bacteria have been shown to 
inject soluble ACs into mammalian cells as a toxin (Sayner 
et al., 2004). Early experimental studies suggested that 
hormone-specific cAMP compartments were generated 
by distinct coupling of receptors and G proteins with 



TABLE 1 

Comparison of computational models that have predicted mechanisms contributing to cAMP compartmentation 



Model Cell type cAMP synthesis cAMP degradation Physical barriers cAMP buffers cAMP export Cell shape 
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subpopulations of transmembrane ACs in specific com- 
partments (Brunton et al., 1981) . Localization of different 
AC isoforms with distinct regulators of cAMP synthesis 
(e.g., (3i-adrenergic receptors, ^adrenergic receptors, G s , 
PKA) at caveolae (Rybin etal., 2000; Ostrom etal., 2001) 
or on specific AKAP complexes (Kapiloff et al, 2014) has 
been shown to be a potential molecular mechanism for 
the formation of specific cAMP microdomains (Fig. 1 B) . 

CNG channels and fluorescent biosensors are fre- 
quendy used to measure localized cAMP synthesis by 
subsets of ACs. Rochais et al. (2004) expressed CNG 
channels in rat ventricular myocytes to compare the 
responses to (3-adrenergic and AC agonists, finding that 
(3-adrenergic receptors activate only ^25% of ACs. This 
measurement quantitatively validated predictions from 
a previous computational model of (3-adrenergic/ cAMP/ 
Ca 2+ signaling (Saucerman et al., 2003). Rochais et al. 
(2006) went on to compare subsarcolemmal cAMP signals 
and changes in L-type Ca 2+ currents, finding receptor- 
specific signatures in response to isoproterenol, glucagon, 
and PGE1. As an alternative approach, the HCN2-camps 
cAMP biosensor showed that local (31-AR-elicited cAMP 
signals propagated globally, whereas cAMP signals gen- 
erated by local (32-AR stimulation remained localized 
(Nikolaev et al., 2006). This approach was elegantly 
extended to also include scanning ion conductance, 
allowing delivery of (3 TAR and (32-AR stimuli specifically 
to t-tubules or outer sarcolemma (Nikolaev et al., 2010) . 
Their data indicated that (32-AR-mediated cAMP syn- 
thesis is normally limited to t-tubules, but these recep- 
tors redistribute to outer sarcolemma in heart failure 
(Nikolaev et al., 2010). The Zaccolo group tethered 
Epac-based cAMP sensors to type I and type II PKA reg- 
ulatory subunits and assessed receptor-specific cAMP 
signals in these compartments (Di Benedetto et al., 
2008). These experiments revealed higher cAMP in the 
type I PKA compartment with (3-AR stimulation, whereas 
glucagon and PGE1 generated higher [cAMP] in the 
type II PKA compartment (Di Benedetto et al., 2008). 

The Harvey group has developed computational 
models that examine the distinct cAMP microdomains 
induced by stimulus-specific local cAMP synthesis. Iancu 
et al. (2007) extended models of (3-adrenergic/ cAMP 
signaling in cardiac myocytes (Saucerman et al., 2003, 
2006) by incorporating caveolar, extra-caveolar mem- 
brane, and cytosolic compartments, each with experi- 
mentally based distributions of receptors ((31-AR and 
muscarinic) and AC isoforms (Iancu et al., 2007). This 
model predicted that baseline caveolar [cAMP] is con- 
siderably below cytosolic [cAMP], preventing constitu- 
tive activation of caveolar PKA. They used this model 
to predict localized cross talk between (31-adrenergic 
and M 2 -muscarinic receptors. Mg-muscarinic receptors 
were predicted to inhibit (3 TAR responses in caveolae 
but enhance cAMP synthesis outside caveolae (Iancu 
et al., 2007) , leading to a quantitative explanation for the 



intriguing acetylcholine rebound effect (Belevych etal., 
2001). Their group subsequendy compared the PKA- 
based cAMP biosensor with a cytosolic Epac-based sensor, 
allowing improved validation of their previous model 
predictions of elevated cytosolic [cAMP] and muscarinic 
stimulation of cytosolic [cAMP] increases (Iancu et al., 
2008) . Later, this compartmental model was extended by 
adding (3 2 -ARs, PKA, and downstream signaling (Heijman 
et al., 2011). Heijman et al. (2011) performed a series of 
virtual knockouts and perturbations, predicting (3 r AR- 
versus pVAR-specific local cAMP signals and distinct 
roles of AC5/6 and AC4/7 regulating caveolar/cyto- 
solic [cAMP] and extracaveolar [cAMP], respectively. 

New experimental techniques for controlling local 
cAMP synthesis have provided additional insights into 
mechanisms of cAMP compartmentation. Saucerman 
et al. (2006) found that receptor-stimulated cAMP syn- 
thesis generated PKA activity gradients between plasma 
membrane and adjacent cytosol, whereas cAMP released 
by local uncaging in the cytosol had a much larger range 
of action. These measurements were quantitatively ex- 
plained using their computational models (Saucerman 
et al., 2006). ACs have also been reengineered to pro- 
vide further control of local cAMP synthesis. To explain 
the paradoxical effects of transmembrane versus soluble 
ACs on endothelial barrier function (Sayner et al., 2004) , 
Sayner et al. (2006) used a cytosolic-targeted AC mu- 
tant, finding that indeed the membrane versus cytosolic 
location of cAMP synthesis determined whether cAMP 
increased or decreased barrier function. 

Development of genetically targetable and bicarbon- 
ate-sensitive soluble AC constructs allowed analysis of 
how cAMP and PKA activity gradients depend on the 
subcellular location of cAMP synthesis (Sample et al., 
2012) . They found that cAMP synthesis at either plasma 
membrane or nucleus generated high local [cAMP] but 
lower distal [cAMP] in HEK293 cells. Further, nuclear 
PKA activity was rapidly activated by nuclear cAMP, in 
contrast to computational models based on the classical 
paradigm of slow nuclear import of PKA catalytic sub- 
unit. Iteration between alternative computational mod- 
els and cAMP biosensor experiments with localized 
cAMP synthesis allowed inference of a resident pool of 
PKA holoenzyme bound to AKAPs and insulated from 
global cAMP by tethered PDE4 (Sample et al., 2012). A 
unique aspect of the modeling approach in this study 
was the iterative rejection and revision of computational 
models together with new experiments to identify a new 
cAMP/PKA microdomain. 

Overall, several computational models have been es- 
tablished to evaluate the impact of localized cAMP syn- 
thesis on generation of [cAMP] gradients. Despite a 
variety of interesting experimental studies showing 
receptor-specific cAMP microdomains, it is important 
to remember that localized cAMP synthesis may not 
fully explain functional differences. Such differences 
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could also be related to the magnitude or kinetics of the 
stimulus, PDE-mediated cAMP degradation, or cAMP- 
independent pathways. Although several of the experi- 
mentally observed receptor-specific cAMP microdomains 
discussed above are not yet accurately predicted, these 
models will serve as a useful framework for future ex- 
perimental and computational studies. Such complexi- 
ties further indicate the need for computational models 
in this area. 

Physical barriers to cAMP diffusion 

Although a wide range of experimental studies have 
clearly demonstrated that cAMP gradients are regu- 
lated by localized cAMP synthesis and degradation, 
even the simplest (and most elegant) computational 
models indicated that this picture may be incomplete. 
Rich et al. (2000) presented a simple analytical model 
of cAMP diffusion from a single AC protein. This analy- 
sis demonstrated that in contrast to the substantial gra- 
dients in Ca 2+ in the vicinity of Ca 2+ channels (Fischmeister 
and Horackova, 1983), biochemically measured rates 
of AC activity are insufficient alone (by a factor of > 100) 
to generate meaningful cAMP gradients (Rich et al., 
2000) . A similar finding was obtained even when ex- 
tending analytical models with localized PDEs (Chen 
et al., 2008). Indeed, only models with physical barri- 
ers to cAMP diffusion between plasma membrane and 
cytosol could predict experimentally measured distinct 
local cAMP gradients induced by forskolin or PGE1 (Rich 
etal., 2000, 2001). 

Several computational models were consistent with 
this concept of restricted cAMP diffusion (Fig. 1 C) . 
Prediction of experimentally measured time delays be- 
tween membrane and cytosolic PKA activity required 
substantially reduced cAMP diffusion rates (by >100x) 
(Saucerman et al., 2006). Simulations of PKA gradients 
generated by local cAMP uncaging predicted that 
cytosolic cAMP diffusion may also be restricted but to 
a lesser degree (Saucerman et al., 2006). Likewise, 
the Iancu et al. (2007, 2008) models predicted that 
low basal [cAMP] in caveolae and compartment-specific 
cross talk between (31-adrenergic and muscarinic signal- 
ing depended critically on restricted cAMP diffusion 
between membrane compartments and cytosol. Similar 
conclusions were obtained in models of distinct cAMP 
microdomains regulated by (31-AR and (32-AR in cardiac 
myocytes (Heijman et al., 2011). The role of restricted 
cAMP diffusion has also been tested in simulations with 
realistic geometries of endothelial cells (Feinstein etal., 
2012). cAMP synthesis and degradation rates were var- 
ied over wide ranges, demonstrating that unless AC/ 
PDE-mediated cAMP synthesis/degradation fluxes were 
highly elevated (by >100x), cAMP gradients required 
restricted cAMP diffusion (Feinstein et al., 2012). 

However, some computational models have indi- 
cated that a reduction in the diffusion coefficient is not 



required for cAMP gradients. The previously described 
models of cAMP gradients in realistic neuron geome- 
tries (Neves et al., 2008) used an unrestricted cAMP 
diffusion coefficient of 200 um 2 /s. cAMP gradients in 
these cells relied on a high PDE activity and additional 
features of cell shape as described below in Regulation 
of cAMP gradients by cell shape. The stochastic models 
of Oliveira et al. (2010) predicted that local cAMP 
synthesis together with broadly distributed and PKA- 
sensitive cAMP degradation were sufficient to generate 
experimentally observed cAMP gradients in HEK293 
cells. The difference in this model may be caused ei- 
ther by the stochastic nature of the model or the assumed 
>500-fold higher cAMP synthesis and degradation rates 
(Vmax PDE, 55-348 uM s -1 ) than had been measured 
biochemically and included in the other models (Rich 
et al., 2000, 2001; Saucerman et al., 2006; Feinstein 
etal., 2012) (see Table 2). 

Although prediction of physical barriers restricting 
cAMP diffusion was certainly a surprising result (Rich 
et al., 2000), evidence of restricted diffusion of other 
small molecules has been reported as well. Fluorescent 
imaging of [Na + ] in ventricular myocytes with local 
Na + /K + pump inhibition identified substantial intra- 
cellular [Na + ] gradients, with an effective diffusion co- 
efficient ^100-fold lower than in aqueous solutions 
(Despa et al., 2004). The distribution of 3-100-nm 
gold nanoparticles injected into permeabilized car- 
diac myocytes indicated substantial physical barriers 
both in the cytosol and subcellular spaces such as 
the dyadic cleft and nucleus (Parfenov et al., 2006). 
Restricted diffusion of nucleotides closely related to 
cAMP has been found as well. Jepihhina et al. (2011) 
measured the kinetics of metabolic respiration in car- 
diac myocytes in response to locally applied ADP, find- 
ing delays in response that could only be explained by 
modeling diffusional obstacles. Correlation micros- 
copy has been applied to fluorescently labeled ATP, 
measuring effective diffusion coefficients of ^4 um 2 /s 
with periodic physical barriers of ^l-um spacing (Illaste 
etal., 2012). 

Overall, most computational models based on bio- 
chemically measured rates of cAMP synthesis and deg- 
radation indicate that physical barriers are necessary to 
predict physiological cAMP gradients. However, the iden- 
tity of these physical barriers is unknown, with possibili- 
ties including endoplasmic/sarcoplasmic membranes 
(Rich et al., 2000), cortical actin cytoskeleton, and myo- 
filament proteins in muscle cells. Although gaps in actin 
cytoskeleton may be considerably larger than cAMP, it has 
been noted that actin networks form gels that consider- 
ably decrease local viscosity (Luby-Phelps et al., 1986). 
A major challenge remains to identify and experimen- 
tally perturb physical barriers to test these model pre- 
dictions without unintentionally manipulating other 
variables such as AC/PDE activity and localization. 
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Restricted diffusion caused by cAMP buffering 
In addition to physical barriers, cAMP-binding proteins 
may bind cAMP to restrict cAMP diffusion (Fig. 1 D). 
PKA is the most likely candidate cAMP buffer, as it binds 
cAMP with higher affinity or is present in larger quantities 
than Epac, PKG, or CNG channels (Rich et al., 2001; 
Poppe et al., 2008) . Early biochemical studies indicated 
that an appreciable fraction of total cAMP is bound to 
PKA, even in unstimulated conditions (Beavo et al., 
1974; Corbin et al., 1977). The model of ^-adrenergic 
signaling developed by Saucerman et al. (2003) ac- 
counted for both free and bound cAMP, predicting 
that PKA buffers the total [cAMP] by shielding it from 
PDE-mediated degradation. cAMP buffering was hy- 
pothesized to stabilize cAMP in the region of maximal 
PKA sensitivity (Saucerman et al., 2003). cAMP buffer- 
ing by PKA was later revisited in the context of localized 
cAMP release, where cAMP buffering was predicted to 
contribute to cytosolic cAMP/PKA activity gradients 
(Saucerman et al., 2006) . cAMP buffering has also been 



predicted to contribute to cAMP gradients in endothe- 
lial cells, although primarily in cases where cAMP buff- 
ering was elevated or coupled with physical barriers 
(Feinstein et al., 2012). Although cAMP buffering has 
been predicted to restrict cAMP diffusion in some cases, it 
should be noted that PKA's buffering capacity for cAMP 
(^1.2 umol/L cytosol; Corbin et al., 1977; Saucerman 
et al., 2003) is much lower than the buffering capacity 
for Ca 2+ . Direct experiments that manipulate cAMP 
buffering while measuring cAMP gradients are needed 
to validate these model predictions and resolve this issue. 

Regulation of cAMP gradients by cell shape 

The first direct measurements of receptor-mediated 
cAMP gradients showed elevated [cAMP] in neuronal 
dendrites, with lower [cAMP] in the cell body (Bacskai 
et al., 1993). Bacskai et al. (1993) suggested that one of 
the simplest explanations for these cAMP gradients 
was the difference in surface-to-volume ratio; dendrites 
with large surface-to-volume ratio would be enriched in 



TABLE 2 

Key numbers for cAMP compartmentation 
Biological parameters Values 

Length scales Adult cardiac myocyte, 120 pm; HEK293, 16 pm; endothelial cell, 40 pm; neuronal dendrite length/ 

width, 100/3 um; cAMP hydrodynamic radius, 1.4 nm 
Concentrations Free cAMP, 0.1-5 |iM; type I PKA, 0.59 uM; type II PKA, 25 nM; Saucerman et al., 2003; CNG channels, 

40 nM; Rich et al., 2001 

Diffusion coefficients for cAMP Aqueous solution, 444 pm 2 s _1 ; Dworkin and Keller, 1977; CHO cell, 487 um 2 s _I (v~ 40 pm s -1 ); Nikolaev 

et al., 2004; Cardiac myocyte, 136 pm 2 s _1 (v= 16 pm s _1 ) ; Nikolaev et al., 2006. Physical barriers between 
membrane and cytosol: 0.5 pm 2 s~ over 1 pm; Rich et al., 2000; 2 pm 2 s _1 over 2.5 pm; Saucerman et al., 
2006 

Rate constants 

AC: kcat-AC, 12-60 s" 1 ; Rich et al., 2000; 8.5 s _1 ; Saucerman et al., 2003 

20 s'\ Feinstein et al., 2012; 

Vmax-AC: cardiac myocyte, 0.1-0.4 pM s -1 ; Saucerman et al., 2003, 2006; endothelial cell, 0.14 pM s" 1 ; 
Feinstein et al., 2012; HEK293, 0.13 pM s -1 ; Rich et al., 2001 
PDE: kcat-PDE3: 3.5 s -1 ; kcat-PDE4: 5 s -1 ; Saucerman et al., 2004; Vmax-PDE: cardiac myocyte, 0.2-0.4 pM s -1 ; 

Saucerman et al., 2003; endothelial cell, 0.3 pM s -1 ; Feinstein et al., 2012; neuron, 3.4 pM s -1 ; Neves et 
al., 2008 

MRP4: k-MRP4: 10"" to 10~ 3 s" 1 ; Cheepala et al., 2013 

Back-of-the-envelope calculations Time for cAMP to diffuse across a cell: 

Cardiac myocyte: 

Aqueous diffusion rate: 4 s 
Cytosolic diffusion rate: 13 s 
HEK293: 

Aqueous diffusion rate: 0.07 s 

Time for cAMP to diffuse from membrane to cytosol in cardiac myocyte: 1.6 s 
Time to degradation of cAMP by PDE: 0.5-50 s 
Distance for cAMP diffusion before degradation by PDE: 

Aqueous diffusion rate: sqrt(2*444*t D ) = 21-210 pm 
Cardiac myocyte cytosol: 12-117 pm 
Restricted diffusion cardiac myocyte = 1-10 pm 
CHO cytosol: 22-220 pm 

Rate constants were obtained from example models as determined from biochemical experiments. Diffusion calculations performed using the solution to 
the 1-D diffusion equation, <x 2 > = 2Dt, where <^> is the mean-squared distance traveled, D is the diffusion coefficient, and t is time (Codling et al., 2008) . 
Using a linear approximation of cAMP degradation by PDE, the time to degradation was calculated by Ueg= [cAMPJ/Vmax — PDE. The time to degradation 
was then used in the 1-D diffusion equation to determine the diffusion distance before degradation. Note that the range of action of cAMP is predicted to 
increase with increasing [cAMP] and decrease with increasing [PDE]. 
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transmembrane AC but contain less cytosolic PDE 
(Fig. 1 E). 15 years later, Neves et al. (2008) performed 
similar measurements with fluorescent protein-based 
cAMP biosensors, again finding elevated [cAMP] in den- 
drites compared with cell body. They built a computa- 
tional model of local cAMP signaling, predicting that 
indeed the compartmental differences in surface-to- 
volume ratio were sufficient to quantitatively predict cAMP 
gradients, without decreasing the cAMP diffusion co- 
efficient (Neves etal., 2008). Chen et al. (2009) further 
developed models of neurons with idealized geometries, 
deriving analytical expressions relating the length and 
width of the dendrite, cAMP degradation rate, and 
[cAMP] gradient magnitudes. The models of Feinstein 
et al. (2012) examined how cAMP gradients may be af- 
fected by physiological variation in endothelial cell shape. 
They found that although surface-to-volume ratios typi- 
cal of endothelial cell geometries were insufficient to 
generate cAMP gradients alone, cell shape did signifi- 
candy enhance the sensitivity to other mechanisms such 
as restricted cAMP diffusion (Feinstein et al., 2012). 

In summary, several computational models predict 
that subcellular geometries with high surface-to-volume 
ratios enhance cAMP gradients, primarily when those 
geometries are at least 5 pm. This indicates that surface- 
to-volume ratio may be insufficient to explain cAMP 
gradients in nanometer-scale compartments such as 
dyadic clefts or caveolae. Experiments that manipulate 
cell shape will be needed to test these model predic- 
tions. But like predictions for physical barriers, it may 
be difficult to perform direct perturbations without also 
changing other key variables. 

cAMP export 

cAMP can also be eliminated from the cell by export via 
multidrug resistance proteins (MRPs) (see Fig. 1 F). Al- 
though MRP-mediated cAMP export rates are very low 
(2 x 10 3s_1 ) relative to overall PDE activities with small 
effects on global [cAMP] (Cheepala et al., 2013), several 
studies have shown that inhibition of MRP can consider- 
ably elevate submembrane [cAMP] (Li et al., 2007; Xie 
et al., 201 1) . However, these studies are complicated be- 
cause pharmacologic MRP inhibitors block PDE activity, 
whereas genetic MRP knockdown can alter PDE expres- 
sion (Xie et al., 2011). Models extended to include 
MRP-mediated cAMP export demonstrated that modu- 
lation of local but not global cAMP by MRP4 required a 
membrane compartment with highly restricted cAMP 
diffusion, low PDE, and high MRP4 (Xie et al., 2011). 

Conclusions and future directions 

The study of cAMP compartmentation has undergone a 
series of paradigm shifts enabled by new technologies, 
from inference by functional measurements to indirect 
biochemical measurements, to direct measurements 
in live cells. Since the advent of cAMP-sensitive CNG 



channels and fluorescent biosensors, considerable ex- 
perimental efforts have been focused on understanding 
the mechanisms underlying cAMP gradients. Of these 
potential mechanisms, local cAMP degradation by PDEs 
has been best characterized because of the availability 
of specific pharmacologic and genetic perturbations. 
Reasonably strong experimental support has also been 
provided for the role of local cAMP synthesis, although 
better targeted tools for measuring and perturbing local 
cAMP signals (e.g., at AKAPs, caveolae) are still needed 
to explain several receptor-specific responses. A variety 
of computational models have been developed to quan- 
titatively explain how previously measured cAMP mi- 
crodomains are regulated by local cAMP synthesis and 
degradation (Table 1). Models have also been used to 
quantitatively explain experimentally observed regulation 
of cAMP microdomains by cell shape and cAMP export. 

In addition to their explanatory power, in several cases 
de novo predictions from computational models have 
been validated in subsequent experiments. Early models 
of (3-adrenergic signaling (Saucerman et al., 2003) ac- 
curately predicted that (31-AR activate only a subpopula- 
tion of ACs (Rochais etal., 2004) , restricting amplification 
of the pathway. Model predictions of cAMP as the rate- 
limiting step in P-adrenergic signaling and predictions 
of time delays between membrane and cytosolic cAMP 
signals (Saucerman et al., 2006) were also validated in 
separate experiments (Leroy et al., 2008). A model of 
(3-adrenergic/M 2 -muscarinic receptor cross talk (Iancu 
et al., 2007) predicted the subsequently validated cy- 
tosolic [cAMP] signals underlying the rebound effect 
caused by transient aceytylcholine (Iancu et al., 2008). 
As a last example, models of nuclear cAMP and PKA 
activity correctly predicted a distinctly regulated nuclear 
PKA-AKAP-PDE4 complex in HEK293 cells (Sample 
etal., 2012). 

Experimental validation of predicted compartmentation 
mechanisms. Additional experimental studies are needed 
to directly validate novel predictions of other cAMP com- 
partmentation mechanisms. Experimental validation of 
the role of cAMP buffering may be feasible by the care- 
ful use of specific cAMP analogues (Poppe et al., 2008). 
In contrast, the identification and specific perturbation 
of physical barriers modulating cAMP compartmenta- 
tion appears considerably more challenging. Further 
experimental validations regarding the role of cell 
shape are also required, including perturbations that 
modulate cell shape dynamically and examination of 
additional cellular structures. For example, PKA activity 
gradients have been observed at the leading edge of mi- 
grating cells (Lim et al. , 2008 ) , which has a much higher 
surface-to-volume ratio than at the rear. 

Modeling functional consequences of cAMP compartmen- 
tation. In addition to inferring biophysical mechanisms, 
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computational models are also valuable for integrating 
distinct biological processes and predicting the multi- 
scale consequences of mechanistic perturbations (Yang 
and Saucerman, 2011). For example, a model coupling 
(3-adrenergic signaling, cellular excitation-contraction 
coupling, and 3-D electrophysiology across the heart wall 
(Saucerman et al., 2004) was able to mechanistically pre- 
dict how human mutations that disrupts a PKA-AKAP 
channel complex lead to T-wave abnormalities subse- 
quendy measured in patients with long QT syndrome 1 
(Darbar et al., 2005). Heijman et al. (2011) also exam- 
ined the impact of localized cAMP/PKA signaling on 
cardiac physiology, predicting in detail how (31- and (32- 
adrenergic stimuli and cross talk with Ca 2+ /calmodulin 
kinase II lead to distinct regulation of action potentials 
and Ca 2+ transients. Neves et al. (2008) predicted and 
then experimentally validated that cAMP gradients in 
neurons help drive gradients in downstream MAPK sig- 
naling related to synaptic plasticity. Greater use of inte- 
grative, multi-scale models would aid the understanding 
of cAMP compartmentation in health and disease but 
also enable in silico design of novel therapeutics based 
on manipulation of local cAMP signals. 

This Perspectives series includes articles by Karpen, 
Kapiloff et al., Rich et al., and Conti et al. 
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